Effect of inter-subsystem couplings on the evolution of composite systems 
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The effect of inter-subsystem coupling on the adiabaticity of composite systems and that of 
its subsystems is investigated. Similar to the adiabatic evolution defined for pure states, non- 
transitional evolution for mixed states is introduced; conditions for the non-transitional evolution 
are derived and discussed. An example that describes two coupled qubits is presented to detail the 
general presentation. The effects due to non-adiabatic evolution on the geometric phase are also 
presented and discussed. 
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PACS numbers: 03.65.Ta, 07.60.Ly 

The study on the adiabaticity of quantal system may 
be traced back to the mid 1980s, when Berry [| con- 
ceived that a quantal system in an eigenstate, adiabati- 
cally transport round a circuit by varying parameters R 
in its Hamiltonian H(R), will acquire a geometric phase 
in addition to the familiar dynamical phase factor. Since 
then geometric phase became an interesting subject and 
has been extensively studied Q and generalized to 
non-adiabatic evolution [§J , mixed states jfj 0, Q , and 
open systems 0. The geometric phase of a composite 
system in particular has attracted a lot of attention for 
its possible applications in quantum information process- 
ing, where the whole set of universal quantum gates are 
achieved based on the Abelian and /or non-Abelian geo- 
metric operations 

mini mill ei ^ 

view of the ge- 
ometric computation, the adiabaticity of the composite 
system is of course an important issue, because it would 
determine how well the system follows the loops. Never- 
theless, thorough studies aimed to address this issue, in 
particular for a composite system with inter-subsystem 
coupling, are still few and certainly not exhaustive |l5j . 

On the other hand, the geometric phase for mixed 
states is a new subject and much remains to be under- 
stood. Uhlmann [f|was the first to address this issue and 
later Sjoqvist et al. formulated it from the viewpoint of 
quantum intcrferometry [7f ■ This formulation is available 
when the system undergoes an unitary evolution. For 
subsystems that compose a composite system with non- 
zero inter-subsystem couplings, however, the evolution of 
each subsystem is not unitary in g eneral. This problem 
was explored in a recent paper |l5( for a very rare situa- 
tion when both the composite system and its subsystems 
evolve adiabatically, but how the subsystems may evolve 
while the composite system transport adiabatically (or 
non-adiabatically) remains an open question. 

In this paper, we will address these issues by investi- 
gating the adiabaticity of a composite system that con- 
sists of two coupled spin-i subsystems or a pair of quan- 
tum bit. We analyze the case where one of the spin-^ 
is driven by a precessing magnetic field, a case of rele- 
vance to Nuclear Magnetic Resonance (NMR) quantum 
computation llfll a s well as to the test of mixed state ge- 
ometric phases • We calculate and analyze the effects 



of spin-spin couplings on the adiabaticity of the compos- 
ite system and its counterpart of the subsystems, four 
different kinds of time evolution are clarified and illus- 
trated, conditions for those evolutions are presented and 
discussed. 

Let a composite system be govern by the Hamiltonian 



H — Hi + H2 + H12 



(1) 



where Hi(i = 1,2) denote the free Hamiltonian of subsys- 
tem i and H12 stands for the interaction between them. 
We suppose that the Hamiltonian H is changed by vary- 
ing parameters R — (X, Y, ...) on which it depends. Then 
the excursion of the system between times t = and 
t = T can be pictured as transport round a closed path 
R(t) in parameter space with Hamiltonian H(R(t)) and 
such that R(T) = R(0). At any instant, the natural basis 
consists of the eigenstates \<fr(R)) of H(R) for R — R(t), 
that satisfy H(R)\(j) n (R)} — E n (R)\(f> n (R)) , with energies 
£ n (R). If H(R) is altered slowly such that 



£n(R) — £m(R) 



« 1, 



(2) 



it follows from the adiabatic theorem that at any instant 
the system will be in an eigenstate of the instantaneous 
Hamiltonian. In particular, if the Hamiltonian is re- 
turned to its original form, the composite system will 
return to its original state, apart from a phase factor. 
Eq. (J2J) is the well known condition for the adiabatic the- 
orem to hold. 

We next develop a generalization for the subsystems, 
going back to the original adiabatic scenario in which 
the system returns to its original state, but now taking 
mixed states into account instead of pure states. To this 
end, we first of all define non-transitional evolution for 
mixed states |18|. this definition is non-trivial in partic- 
ular for subsystems that have no effective Hamiltonian 
available for it[l9|. Let a state p(t) of the subsystem 
(say, subsystem 1) be written in the diagonal form of 
p(t) = 52iPi(t)\Ei(t))(Ei(t)\, p{t) depends on time via 
R(t) and we would write the time-dependence of R(t) 
explicitly. It is clear that Pi(t) gives the probability of 
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the subsystem in state \Ei(t)). This form of writing is 
called the spectral representation, while Pi (i) denotes the 
eigenvalues and \Ei(t)) the corresponding eigenvectors 
of p(t). One special case is that Pi{t) (for any i) are 
time- independent, this is a rare situation which implies 
no transitions among the eigenstates of p(t) when the 
composite system experiences transport along the pa- 
rameter loops. The subsystem in this state with time- 
independent coefficients pi (t) , i.e., Pi(t) = Pi(0) indepen- 
dent of the varying parameters, is defined to undergo 
non-transitional evolution, and the corresponding eigen- 
states \Ei(t)) will be called non-transitional eigenstates. 
Obviously the non-transitional evolution would return to 
the adiabatic evolution when the states \Ei(t)) are the 
eigenstates of the subsystem's Hamiltonian(if any avail- 
able). Moreover this definition is meaningful even if there 
is no Hamiltonian available for the subsystems, a gen- 
eral situation for coupled multi-particle systems. Thus 
the definition could find broad use instead of the adia- 
batic evolution for pure states in composite systems. Now 
we will drive a condition for the subsystem to undergo 
this non-transitional evolution. For the composite sys- 
tem govern by Hamiltonian Eq.©, a state may be 
decomposed into Schmidt form 



J2VpW)^ iI ° Hu(t ' )dt ' mt))M(t)h, (3) 



with I 2 (e l (t)\H\e J (t)) 2 \E J (t)) 1 and 

l e j(*)>2j-Sj(*)>i = |ej(*)>2 ® \Ej{t))x. The reduced 
density matrix for the subsystem 1 follows straightfor- 
wardly from Eq.© that p x {t) = Tr 2 (|^(i)) (V>(*)|) = 
Y /l Pi(t)\E l (t))(E t (t)\. To find the condition of the 
non-transitional evolution is now equivalent to finding 
conditions for p%{t) in Eq.© to be time-independent. In 
units with K = 1, \ip{t)) satisfies 



i^\m) = H\m), 



(4) 



where here and hereafter time-dependence are under- 
stood where not written explicitly and \Ei(t)) (\ei(t))) 
denotes states for subsystem 1 (2) where suffix omitted. 
We would like to note that \Ei(t))\ei{t)) are not the in- 
stantaneous eigenstates of H in general, so under the ac- 
tion of H transitions among those states would occur. As 
you will see, the condition for the non-transitional evolu- 
tion would be equivalent to negligible ratios of the tran- 
sition amplitude to the respective energy spacing. The 
derivative equation for pi{t) follows from Eq.© that, 

iy/Pj + VPj( E j\ E j) + \/Pj< e i|ej) 
- ^2 Vpk~exp{-i / (Hkk - Hjj)dt }H ]k = 0. (5) 

The simplest approximation is to neglect the off-diagonal 
elements on the grounds that 



having this approximation, Eq. © yields 

y / ^= v /^e^+^), (7) 

where 7y = i ^{E^^E^dT, and -y 2 j = 

iJo( e j( T )\w\ e j( T )) dT - We may rewrite |£y(r)) (Mt))) 
to be e^\Ej(T)) (e 1 ^ 2 |ej(r)» such that pj(t) = Pj (0). 
Usually, \Ej(t))\ej(t)) are not the instantaneous eigen- 
states of the Hamiltonian H, so Hjk represent the 
transition amplitude between states \Ej(t))\ej(t)) and 
\Ek(t))\ek(t)) . Condition Eq.© indicates that for non- 
transitional evolutions, transitions induced by H among 
|-^t(*))| e i(*)) should be small with respect to the energy 
spacing between the two. Conditions © and © together 
imply that there might be four different kinds of evolution 
for the composite system, as follows, (a) Adiabatic evolu- 
tion: The composite system undergoes an adiabatic evo- 
lution while its subsystems follow non-transitional evolu- 
tions, this means that the composite system would follow 
one of its instantaneous eigenstates while its subsystem 
evolve along the non-transitional eigenstates. (b) Quasi- 
adiabatic evolution 1: The composite system under- 
goes an adiabatic evolution while its subsystems do not; 
(c) Quasi- adiabatic evolution 2: The composite system 
evolve non-adiabatically while its subsystems follow non- 
transitional evolutions. (d)N on- adiabatic evolution: The 
composite system evolve non-adiabatically, while its sub- 
systems undergo out of non-transitional evolutions. 

It is worthwhile to mention that the condi- 
tion/criterion for non-transitional evolution given by 
Eq.© is also valid for the system with a mixed state 
under unitary evolutions, this can be understood as fol- 
lows. Write the state (density matrix) of the composite 
system in a spectral representation 



p(t) = }^\ a \^ a (t))(ip a (t)\, (8) 

a 

\ip a (t)} would obey the Schrodingcr equation and has the 



same form of decomposition as Eq.© but with yp^(t) 

instead of y/pj(t). The von Neumann equation ift^f = 
[H, p] then yields the same equation as in Eq. © but with 



^fp~j = y X)a-Pj ^a- J2P*j^a represents the population of 
the subsystem 1 in state \Ej(t)}, non-transitional evolu- 
tion that requires pj constant would result in the same 
condition as in Eq.©, indeed for a closed system the 
transitions among the eigenstates of the subsystem's den- 
sity matrix are only induced by the system Hamiltonian, 
for open systems, however, this is not the case, since the 
bath which couples to the whole system may lead to pop- 
ulation transfer among those states. To this respect, in 
addition to the condition Eq.©, criteria 

,(E fc |(e fc |r a |e 3 )|E/ 



« 1 



(9) 



«1, 



(6) 



is required to be satisfied in order to ensure the non- 
transitional evolutions for open systems, where T a rep- 
resent the agents (operators) of the whole system coupled 
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to a bath. This is the reason why the non-transitional 
condition should be different from each other for an open 
system and for a closed system, and this is also an expla- 
nation for closed systems, regardless of in a pure state or 
a mixed state, hold the same criteria for non-transitional 
evolutions. 

As an example, we consider two qubits Sk(k = 1, 2) as 
represented by a pair of spin-i particles, coupled through 
a uniaxial exchange interaction in the z-direction. One 
of the qubits (say, qubit 1) is driven by a time-dependent 
magnetic field B(t) — B h(t) with the unit vector n = 
(sin 9 cos 4>, sin 9 sin e/), cos 9), the Hamiltonian of this sys- 
tem reads (h = 1) [2(j 



H (t) = USf <g> S$ + (iB(t) • Si 



(10) 



with the exchange interaction constant J and the gyro- 
magnetic ratio /i. This Hamiltonian is of relevance to 
NMR experiment where Carbon- 13 labelled chloroform 
in de acetone may be used as the sample. The single 13 C 
nucleus and the 1 H nucleus play the role of the two spin- 
i particles; the spin-spin coupling constant in this case 
is 4J ~ (27r)214.5Hz. The individual addressing can be 
realized in NMR by using spins of nuclei of different iso- 
topes, such as those of different species of atoms, as the 
subsystems, these spins usually have precession frequen- 
cies that differ from the other by many MHz, a resonant 
magnetic field for one spin then has little effects on the 
others. For more detail, we refer the reader to Ref . plj . 

The instantaneous eigenstates and the corresponding 
eigenvalues (in units of ^/j,Bq) can be written as 



M*)) 
M*)> 



i 



i 



3,4 



(sin0e-^| iT> + (s + cos0 + £i, a )| ft)), 

(sin0e-^| U) + (cos + £ 3 ,4 -5)1 U», 
I 



(11) 



and 



units of fiB /2) ui << 



I- 4 - 

1 sin ( 



Obviously, in this limit 



£1,2 = ±v/( 3 2 + l) + 2 5 cos0, 
£3,4 = ±V(9 2 + l)-2 5 cos0, 



(12) 



where g = denotes the rescaled exchange inter- 

action constant, Mj the renormalization constant and 
T4) = I T)i ® I D2 and the others likely. For a 
simplest case where the exchange interaction g = 0, 
the eigenvalues are reduced to £± = ±1, with cor- 
responding instantaneous eigenstates = (cos | 
I T) +sinfe-^) ® I T) (or ®| |)), |0_(t)) = (- sin | 
I T> + cosfe"^) ® I t) (or ®| 4)). Suppose the ex- 
ternal magnetic fields precess with time-independent az- 
imuthal angles 9 and a constant precessing frequency u, 
i.e., 4> — tot, the adiabatic evolution for the composite 
system requires |(0+| </>_)/(£+ — £-)\ << 1, that is (in 



the two spin-^ particles remain uncoupled and the adi- 
abatic condition is exactly the one for the driven parti- 
cles 1. From the aspect of non-transitional evolution, no 
constraints on to could be made, because of there is no 
coupling between the two qubits and each qubit would re- 
main in pure states if the initial states are pure. So, for 
a composite system without inter-subsystem couplings, 
the evolutions would fall in regime (a) or (c), i.e., the 
composite system might undergo an adiabatic or non- 
adiabatic evolution, while its subsystems evolve along the 
non-transitional states certainly. 



Now we turn to study the case with inter-subsystem 
couplings. For the composite system, to make the adi- 
abatic theorem valid, it should be satisfied that (i,j — 



J 



Tij = 0, others. (13) 



This condition follows straightforwardly from Eq. (J2J) by 
assuming the azimuthal angle 9 time-independent and 
4> = tut. Clearly, the eigenenergies £i and the renor- 
malization constant Mi are independent of u, so in- 



crease linearly with uj, i.e., slowly precessing magnetic 
fields would benefit the adiabatic evolution. The de- 
pendence of = 1,2 or 3,4) on g and 9 was 
illustrated in figure 1. A common feature of these fig- 
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34 




*ft&|*(t)> - H\9(t)) yields (with fi jfc (t) = J*S k (r) - 

= ->Je^W(^(t)|^(t))c fc (t), (16) 

the Berry's phase is just a follow up of this equation by 
ignoring its right-hand side, further consideration would 
treat the terms in the right-hand side as perturbations. 
Up to the first order in I l&Ml^IQ) I the usual perturba- 
tion theory give( assuming c„(i = 0) = 1) 



FIG. 1: (Color online) Selected results for Tij as a function of 
the azimuthal angle 8 [Arc] and the rescaled coupling constant 
g (in units of fiBo/2). u = 10(in units of fiBo/2) was chosen 
for this plot. Figures (a),(b) are for different % and j; (a) Ti2, 

(b) r 34 . 

ures is that Ty — ► with the rescaled coupling constant 
g — ► oo. This limit corresponds to the case when the 
second term in the Hamiltonian Ea. (|10|l can be ignored. 
Physically, the inter-subsystem coupling increase the en- 
ergy spacing between any two instantaneous eigenvalues, 
this makes the population transfer between the respec- 
tive instantaneous eigenstates more difficult, and equally 
the adiabatic evolutions easier. r 12 , and 1^34 behave as 
a non-monotonic function of g as figures l-(a), and (b) 
show, they increase to a maximum value for a disaster 
coupling and then towards to zero for a sufficiently large 
coupling. Two singular points can be found in figure 1, 
for Ti2 the point is located at g = 1 and 9 = ir, while for 
rW, 5 = 1 and 8 = 0. This singularity results form the 
degeneracy of the corresponding eigenstates, and hence 
the adiabatic condition is difficult to meet at these points. 

These transitions induced by non-adiabaticity among 
the instantaneous eigenstates of the Hamiltonian would 
affect the geometric phase of the composite system, the 
geometric phase in this case may be calculated by remov- 
ing the accumulation of these dynamical phases from the 
total phase, i.e., 

4> g =arg{*(p)\9(T))+i C dt(* (t)\* (t)) , (14) 

Jo 

where |^(t)) represents the state of the composite system 
at instance t, and can be written as 

|*(*)) = Ec^e-^^l^)), (15) 

with \4>j(t)) the instantaneous eigenvector and £j the 
corresponding eigenvalue. The Schrodinger equation 



c n (t) ~ e^«, 

^ e^(^^( t )l^( t )) |i/n 
£j - tk 

with J p (t) = i J (<f>p(T)\^> p (r))dT, the Berry phase per- 
tained to instantaneous eigenstate \<f> p (t)). Ea. p4(l and 
ff 1 71) together yield the geometric phase 

4> g = arg^W + ]T c m (T)]. (18) 

To get this expression, all terms equal or smaller than 
T^ nn have been ignored. It is clear that the main cor- 
rection caused by the non-adiabatic evolution to the ge- 
ometric phase come from terms J2 m: £ n Cm(T), that may 
be expressed up to first order in T mn as 

4>g = Jn(T) + ^ T mn[^n m (T) +7 m (T) 

+ arg(4> n (T)\4> m (T))}, (19) 

with r„ m — > the geometric phases approach the Berry 
phases J n (T) as expected. We may divide the correc- 
tion due to the non-adiabatic evolution to the geometric 
phase into two kinds, the first is the population trans- 
fer among the instantaneous eigenstates, this contribu- 
tion appears in the correction as 53m/n T mn (fl nm (T) + 
arg(</) n (T)\(/) m (T))), while the second that exhibit the 
geometric feature of the transited eigenstates scales as 
Z)m#« r ™™7m(r). We would like to note that only T 12 
and ]?34 are not zero in the example Eq. i|10|) . but the 
representation Eq. QlSjl is quite general for quantum sys- 
tems. 

It is not difficult to show from Eq. l|ll|) that the trans- 
port of the subsystems is always non-transitional in this 
situation. For example, suppose the composite sys- 
tem undergo an adiabatic evolution in the instantaneous 
eigenstate \<fr n (t))(n — 1,2), the reduced density matrix 
of the subsystem 1 reads 
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1 ( (<? + cos(9 + £„) 2 
M n \sva6(g + £ n +cos6)e-^ 

P+\p+)(p+\ + P-\P-)(P-\ 



oi<t> 



sin 2 9 



( Pll Pl2 
\ P*12 P22 



(20) 



with /?+.- = 0, 1. Clearly, p± are time- independent. This 
point will be changed when the spin-spin coupling takes 
the form of (Ja^a^ +h.c), the instantaneous eigenstates 
in this situation are entangled states of the two particles, 
the diagonal elements of the reduced density matrix p\ 
is ^-dependent as Ref . shown, and the slowly varying 
9 can make the composite system an adiabatic evolution 
meanwhile make the subsystems out of non-transitional 
evolution. 

Some remarks on the non-transitional evolutions are 
now in order. For subsystems with an available Hamil- 
tonian, the concept of non-transitional evolution covers 
the concept of adiabatic evolution, this can be under- 
stood as follows. Write an initial state (generally mixed) 
of the subsystem in the natural basis (the instantaneous 
eigenstates of the subsystem's Hamiltonian) pi(0) = 
Pij\^i(°))(^j(°)\dh3 = l>-)4), adiabatic evolu- 



tion yields Pl (t) = E^P^WK$;WI,(M = l,-,4). 
Diagonalizing p(t), we can get non-negative and time- 
independent eigenvalues since = p*,, this means that 
the usual adiabatic evolution must fall in the regime of 
the non-transitional evolution, but the inverse could not 
be proven correct. 

In conclusion, the evolution of composite systems and 
its subsystems has been studied. By the definition of 
non-transitional evolution, four different kinds of evolu- 
tion were identified and illustrated via the coupled two- 
qubit system. The non-transitional evolution would find 
its use in formulating evolution of composite systems, 
in particular for subsystems that have no Hamiltonian 
available. 
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